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ABSTRACT 

We present searches for molecular gas (CO, OH, CS, and NH 3 ) in six GHz 
Peaked Spectrum (GPS) radio sources. We do not detect gas in any source 
and place upper limits on the mass of molecular gas which are generally in the 
range ~ 10 9 to a few xlO 10 M . These limits are consistent with the following 
interpretations: (1) GPS sources do not require very dense gas in their hosts, & 
(2) The GPS sources are unlikely to be confined by dense gas and will evolve to 
become larger radio sources. 

Subject headings: galaxies: active - galaxies: ISM - radio lines: galaxies 



1. Introduction 

The GHz Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS) radio sources 
make up significant fractions of the extragalactic bright (cm wavelength selected) radio source 
population (~ 10% and ~ 20%) respectively, but are not well understood (e.g., O'Dea 1998). 
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They are powerful but compact radio sources whose spectra are generally simple and convex 
with peaks near 1 GHz and 100 MHz respectively. The GPS sources are entirely contained 
within the extent of the nuclear narrow line region (< 1 kpc, NLR) while the CSS sources 
are contained entirely within the host galaxy (< 15 kpc). 

GPS and CSS sources are important because (1) they probe the NLR and interstellar 
medium (ISM) of the host galaxy and (2) they may be the younger stages of powerful large- 
scale radio sources - giving us insight into radio source genesis and evolution. 

There are two main hypotheses for the GPS and CSS sources. They could be the young 
progenitor of the large scale powerful double sources (e.g., Carvalho 1985; Hodges & Mutel 
1987; Begelman 1996; Fanti et al. 1995; Readhead et al. 1996b; O'Dea 1998; Snellen et al. 
2000; Alexander 2000). In this model they propagate relatively quickly through the ISM 
of the parent galaxy with advance speeds of a few percent of the speed of light. (Observed 
proper motions tend to be a bit higher - in the range 0.1 - 0.2, e.g., Polatidis & Conway 2003, 
though the detections may be biased towards objects with the highest velocities). In order 
to allow these high velocities, the ISM cannot be very dense and the total cold gas content 
can be no more than about 10 10 M . These sources must also undergo strong luminosity 
evolution as they evolve, dimming by 1-2 orders of magnitude in radio flux density. 

In the second hypothesis, these sources are older and propagate much more slowly 
though a dense ISM acquired via cannibalism (O'Dea, Baum, & Stanghellini 1991; De Young 
1993; Carvalho 1994, 1998). These frustrated sources interact strongly with their dense 
ambient medium driving a shock which ionizes the gas and produces two effects - (1) free 
free absorption which is responsible for the turnover in the radio spectrum and (2) optical 
emission lines (Bicknell, Dopita, and O'Dea 1997). In this second model there is very little 
luminosity evolution. Bicknell et al. predict that the host galaxies will be relatively gas-rich 
with total masses of cold gas in the range 10 10 — 10 11 M . We note that this is similar to the 
cold gas content of the Ultra Luminous Infrared Galaxies - ULIRG (e.g., Sanders, Scoville, 
& Soifer 1991). So far molecular gas (as traced by CO) has been detected in only one GPS 
source (1345+125 - Mirabel, Sanders, & Kazes 1989; Evans et al. 1999b) with estimated 
mass of 6 x 10 10 M . No CO was detected in the GPS source 1934-638 at a limit of 5 x 10 10 
M (O'Dea et al. 1994b). 

Searches for the 21 cm line have produced a 50% detection rate in GPS and CSS sources 
while the detection rates for large sources are only ~ 10% (Vermeulen et al. 2003; Pihlstrom, 
Conway & Vermeulen 2003; van Gorkom et al. 1989). This indicates that clouds of atomic 
hydrogen are very common in the environments of GPS and CSS sources. Recent Hubble 
Space Telescope imaging and spectroscopy has shown high surface brightness emission line 
gas is aligned with the radio axis in CSS radio galaxies (De Vries et al. 1997,1999, Axon et 
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al. 2000; O'Dea et al. 2002). 

Thus, the two models for GPS/CSS sources predict a substantial and testable differ- 
ence in the cold gas content of the host galaxies. And there is some existing evidence that 
GPS/CSS sources contain dense gas in their host galaxies. We have undertaken two com- 
plementary searches for molecular gas in GPS sources. First, we obtained VLA observations 
to search for several molecular species (NH 3 , CS and OH) in four GPS sources. Second, we 
obtained IRAM 30m observations of three (relatively) low-redshift GPS sources to detect 
or set limits on the cold gas content. These three sources are the lowest redshift sources 
from volume limited subsets of complete samples of GPS (Stanghellini et al. 1998) and CSS 
(Fanti et al. 1990) sources. One object, 0428+205 is observed in both the VLA and IRAM 
searches. 

In this paper we present searches for molecular gas with the VLA and IRAM in 6 
GPS sources. In principle, both emission and absorption searches are possible. Given the 
redshfits of the sources and the cm wavelength flux densities, absorption and emission limits 
are most sensitive for the VLA and IRAM data, respectively. We estimate constraints on 
the molecular gas content and discuss the implications for models of these sources. 

2. Observations and Reduction 

2.1. VLA Observations 

We searched for absorption in a transition of either ammonia, CS or OH (whichever 
fell in a VLA band) in four GPS sources (Table 1). We observed with the VLA (Napier 
et al. 1983) on January 22, 1993 in the A configuration in spectral line mode (Mode 1A) 
using on-line hanning smoothing, 32 channels, channel spacing of 390.625 kHz, and 12.5 
MHz total bandwidth for each observation. For the 15 GHz observations of 0237-233 we 
obtained additional observations with the central frequency offset by ±10 MHz to cover a 
total bandwidth of about 33 MHz. Observational parameters are given in Table 2. Bandpass 
calibration was performed using observations of the closest of either 3C84 or 3C454.3. Flux 
density calibration was carried out using observations of 3C48. The data were reduced in 
AIPS following standard procedures. 

A few very noisy channels at both ends of the spectra were deleted. Since the sources 
are compact and unresolved at the VLA resolution we used the task POSSM to average the 
data for all the antennas to obtain the integrated spectrum. We subtracted a least squares 
linear fit to the continuum. 
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2.2. IRAM Observations 

We used the IRAM 30m millimeter-wave telescope, located on Pico Veleta, Spain, to 
search for CO emission in the GPS sources 0116+319, 0428+205, and 0941-080. The 
observations took place over 11-13 July 1998 during daylight hours. Receivers were tuned 
to the redshifted CO 1 — > transition (115 GHz rest frequency) and, for 0116+319 only, 
separate receivers were tuned simultaneously to the CO 2 — > 1 transition (230 GHz rest 
frequency). The beamsize at these transitions is 21" and 11", respectively. The telescope 
pointing and focus were calibrated against scans of Jupiter, Mars, and the BL Lac object 
0235+164. To stabilize the spectral baselines and perform initial sky subtraction, the targets 
were observed with a wobbling secondary. The secondary throw angles ranged from 60" to 
150", and the wobble frequency was 0.25 Hz. 

Each transition was observed in two backends, an autocorrelator with 1.25 MHz channel 
separation over a 600 MHz bandwidth, and a filterbank with 1 MHz channel separation over 
a 512 MHz bandwidth. Spectral baselines were subtracted using low-order polynomial fits to 
the raw spectra. The baseline-subtracted spectra for each source were then averaged using 
statistical (1/RMS 2 ) weighting. The averaged spectra were finally converted to mJy from 
Ta* using nominal sensitivity curves provided by IRAM. The basic properties of the reduced 
spectra are listed in Table 3. 



3. Results 
3.1. VLA Results 

We did not detect any significant absorption in these sources. VLA Results are given in 
Table 4. We obtain upper limits to the absorption optical depth of typically a few percent 
with a total range of 1-10%. 

The column density for OHA1667 absorption is given by 

N(OH) ~ 2.2 x 10 14 T cx rAV cm" 2 (1) 

where T ex is the excitation temperature, r is the optical depth, and AV is the FWHM of the 
line in km s" 1 (for 0428+205 and 2352+495 we adopt the values AV = 297 km s" 1 and AV = 
82 km s -1 , respectively, matching the widths of the broadest HI absorption features found 
by Vermeulen et al. 2003). Following O'Dea & Baum (1987) we adopt a fiducial excitation 
temperature T cx = 10 K (e.g., Dickey, Crovisier, & Kazes 1981; Turner 1985). We convert 
from OH to H 2 column assuming a relative abundance ratio 10~ 7 (e.g., Guelin 1985; Irvine 
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et al. 1985). The mass of molecular gas is given by M(mol) = l.3GirR 2 m H 2N(H2) where 
m#2 is the mass of a hydrogen molecule, R is the radius of the region considered (assuming 
for simplicity a plane parallel geometry) and the factor of 1.36 includes the contribution of 
He to the total molecular mass at solar abundance (e.g., Sanders et al. 1991). 

We applied a similar approach to evaluate the detection limits on the CS and NH 3 
absorption. The column density of CS can be estimated by 

iV(CS) ~ 3.6 x 10 12 [1 - exp(-2.35/T cx )]TAV cm" 2 (2) 

(e.g., Turner et al. 1973; Gardner & Whiteoak 1978). The relative abundance ratio CS / H 2 
in dense Galactic cloud cores is ~ 1CT 9 (e.g., Snell, Langer, & Frerking 1982). As for OH, 
we also adopt T cx = 10 K. 

Only 0237—233 was searched for CS absorption (or emission). The relatively narrow 
velocity coverage of the CS observations limits the detection of absorption lines to line widths 
AV < half of a single bandpass. For the purposes of evaluating molecular mass limits, we 
therefore adopted AV = 0.5 x the velocity range spanned by an individual bandpass, or 
~ 120 km s -1 . In principle, a large amount of molecular gas could "hide" in a larger line-of- 
sight velocity dispersion, but the range of velocities that absorption selects is restricted by 
the compact and narrow background continuum source. The velocity widths of absorption 
lines detected in active spirals and ellipticals tend to be < 150 km s _1 (e.g., Vermeulen et 
al. 2003; Gallimore et al. 1999; van Gorkom et al. 1989; Dickey 1986; Kazes & Dickey 
1985), although very broad absorption lines have been detected (albeit rarely) in systems 
with larger background radio sources (e.g., PKS 2322-123, FWHM ~ 735 km s" 1 - O'Dea, 
Baum, & Gallimore 1994; Taylor et al. 1999). It seems therefore unlikely that large molecular 
columns are suppressed by line-of-sight velocity dispersions (greatly exceeding 150 km s _1 ) 
in this particular compact radio source. 

Assuming LTE, the column density of NH 3 is given by 

iV(NH 3 ) ~ 2.8 x lO^/TiVnAV cm" 2 (3) 

(e.g., Batrla, Walmsley, & Wilson 1984; Herrnstein 2003) where r n refers to the peak optical 
depth of the (J, K) = (1,1) rotation inversion transition, and fn is the fraction of NH 3 
molecules in the (1,1) state. Typical excitation temperatures in Galactic cloud cores are 
T ex ~ 40 K (e.g., Morris et al. 1973; Barrett, Ho, & Myers 1977; Ho et al. 1977), at which 
f u ~ 0.3 (Herrnstein 2003). We adopt the relative abundance ratio NH 3 / H 2 = 10~ 8 (e.g., 
Morris et al. 1973; Turner 1995). Only 0404+768 was searched for NH 3 absorption (or 
emission), and for the purposes of evaluating the molecular mass limit, we assume AV = 
107 km s _1 , based on the broadest HI absorption line detected in this source (Vermeulen et 
al. 2003). 
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3.2. IRAM Results 

The results of the IRAM observations are summarized in Table 5. We do not detect 
CO in any of the 3 GPS sources searched. We obtain an upper limit to the integrated flux 
density of the CO line Wco an d convert this to an upper limit on the molecular gas content 
using the relation given by Sanders, Scoville & Soifer (1991) 

M(H 2 ) = 1.18 x 10 4 ^co i->oAl m M (4) 

where Wco is i n units of Jy km s _1 and the luminosity distance D\ um is in units of Mpc. 
The conversion to molecular hydrogen is subject to systematic uncertainties which have been 
thoroughly discussed elsewhere (e.g., Bloemen et al. 1986; Scoville & Soloman 1987; Young 
& Scoville 1991; Sanders et o/.1991; Maloney & Black 1988; Sage & Isbell 1991). The limits 
derived from the CO 2 — > 1 measurements assume CO 2 — ► 1 / CO 1 — > = 0.6 (Lazareff et 
al. 1989). 

4. IMPLICATIONS OF THE LACK OF DENSE MOLECULAR GAS 

We did not detect molecular gas in any of the 6 GPS sources we searched. The estimated 
upper limits to the masses of molecular gas in the range few xlO 8 to a few xlO 10 M Q 
(Table 5), with the weakest limit of < 10 11 M based on the non-detection of NH 3 absorption 
towards 0404+768. These are less than the amounts detected via CO in the most gas rich 
systems such as ULIRGs (e.g., Sanders et al. 1991). However, most of these GPS sources 
could still possess gas masses consistent with the lower end of the range observed in luminous 
Infrared Galaxies. 

The detection of 21 cm H I absorption toward roughly half of the known GPS sources, 
but the absence of molecular line absorption, suggests two possibilities: either (1) the cir- 
cumnuclear gas is predominantly atomic, or (2) atomic gas is easier to detect than molecular 
gas. For absorption experiments against a very bright continuum source, only the line optical 
depth determines the absorption signal strength. Assuming LTE and large densities, we can 
compare the column densities detectable by the two methods. 

For atomic hydrogen 21 cm line, we can estimate the column density by 

N H = 2 x 10 20 cm" 2 (r/0.1)(A\//10rmA;ms" 1 )(T/100K) 

where Nh is the neutral atomic hydrogen column density, r the optical depth, Ay the 
linewidth, and T the gas temperature. 
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For the CO 1-0 molecular line, the column density is 

N H2 = 5 x 10 21 cm~ 2 (r/0.1)(A\//10rmA;ms~ 1 )(T/100is:) 

Here, we assume that the CO/H 2 abundance ratio is 10 -5 . Results for the other 
molecules in our survey will give similar results. Thus, these estimates show that absorption 
atomic hydrogen is indeed easier to detect than molecular line absorption. 

The few powerful extended radio galaxies which have been detected in molecular gas 
tend to have gas masses in the range 10 9 - 10 10 M (e.g., Israel et al. 1991; Mazzarella et al. 
1993; O'Dea et al. 1994a; Evans et al. 1999a; Lim et al. 2000; Das et al. 2004, in preparation). 
This suggests that the GPS and larger radio galaxies may have similar molecular gas content. 

These results suggest the following implications for GPS sources. 

• GPS sources do not require extremely dense environments. 

• The lack of very large masses of dense gas is consistent with the hypothesis that the 
majority of GPS sources are generally not frustrated and will likely expand to become 
CSS sources. This is consistent with the proper motions of ~ 0.1 — 0.2c observed in 
about ten GPS sources so far (e.g., Owsianik & Conway 1998; Owsianik, Conway, & 
Polatidis 1998; Tschager et al. 2000; see compilation by Polatidis & Conway 2003). 
However, these results do not require that all GPS sources evolve to become large 
classical sources - though evolution is favored by the existing data (e.g, O'Dea 1998). 
Alternate models are still possible for part of the GPS source population - e.g., some 
GPS sources may be intrinsically short lived (Readhead 1995). 
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Table 1. Source List 



Name 


Other Name 


ID 


z 


Dl um 


Transition 


Rest Frequency 










Gpc 




GHz 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


0116+319 


4C31.04 


G 


0.0598 


0.26 


CO 1-0 


115.271 












CO 2-1 


230.542 


0237-233 


OD263 


Q 


2.223 


17.9 


CS 1-0 


48.991 


0404+768 


4C76.03 


G 


0.59846 


3.51 


NH3(1,1) 


23.694 


0428+205 


OF247 


G 


0.219 


1.07 


OH 


1.667 












CO 1-0 


115.271 












CO 2-1 


230.542 


0941-080 




G 


0.2280 


1.12 


CO 1-0 


115.271 












CO 2-1 


230.542 


2352+495 


OZ488 


G 


0.23831 


1.18 


OH 


1.667 



Note. - (1) B1950 IAU Source Name. (2) Other Name. (3) ID, G=Galaxy, 
i.e., narrow lines only, and Q=quasar, contains broad emission lines. (4) Redshift. 
0404+768 and 2352+495 taken from Lawrence et al. (1996). 0428+205 from O'Dea 
et al. in preparation. (5) The luminosity distance estimated using Ned Wright's 
cosmology calculator applet assuming the current A cosmology, i.e., H Q = 71 km s -1 
Mpc -1 , and Vl M = 0.27. (6) The molecular transition observed. (7) Rest frequency 
of transition. 



Table 2. VLA Observation Parameters 



Name 


Integration Time 


Central Frequency 


Velocity Resolution 


Velocity Coverage 




Min 


GHz 


km s 


km s 


0237-233 


39 


15.2004 


7.7 


246 


0237-233 


20 


15.2104 


7.7 


246 


0237-233 


24 


15.2204 


7.7 


246 


0237-233 Total 


83 


15.2104 


7.7 


650 


0404+768 


43 


14.82296 


7.9 


253 


0428+205 


44 


1.36781 


85.5 


2740 


2352+495 


37 


1.34790 


86.7 


2780 



Note. — (1) Source Name. (2) Integration time in minutes. (3) Central frequency in GHz. (4) Velocity 
spacing of a 390 kHz channel. (5) Total velocity coverage in km s _1 . 
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Table 3. IRAM Observation Parameters 



Source 


Backend 


Transition 


Int. Time 


T A * rms 


RMS 


Channel width 










min 


mK 


mJy 


km s _1 


0116+319 


1 MHz Filterbank 


CO 1 - 


-> 


309 


2.0 


12.6 


2.9 






CO 2 - 


-> 1 


309 


5.5 


52.8 


1.5 




Autocorrelator 


CO 1 - 


■+ 


309 


3.4 


21.4 


3.7 






CO 2 - 


-> 1 


309 


4.0 


57.6 


1.8 


0428+205 


1 MHz Filterbank 


CO 1 - 


-> 


720 


1.6 


9.6 


3.9 




Autocorrelator 


CO 1 - 


-> 


720 


1.5 


9.0 


4.8 


0941-080 


1 MHz Filterbank 


CO 1 - 


-> 


610 


2.2 


13.2 


3.9 




Autocorrelator 


CO 1 - 


-> 


610 


2.0 


12.0 


4.9 



Note. - 



Table 4. VLA Results 



AT 

Name 


Flux Density 


rms noise 


j_ * 1 1 j_ 1 

optical depth 




Jv 


mJv 




0237-233 


1.46 


2.3 


0.0047 


0237-233 


1.46 


5.0 


0.0103 


0237-233 


1.46 


3.4 


0.0069 


404+768 


1.45 


2.7 


0.0056 


0428+205 


3.84 


1.5 


0.0012 


2352+495 


2.71 


2.2 


0.0024 



Note. - (1) Source Name. (2) Flux density in Jy. (3) 
channel-to-channel rms noise in mJy For 0237-233 the noise 
is given for each of the three observations in order of in- 
creasing central frequency. (4) 3cr upper limit on the optical 
depth of the line. 
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Table 5. Constraints on Molecular Gas 



Name 


Transition 


M(H 2 ) 






M© 


(1) 


(2) 


(3) 


0116+319 


CO 1-0 


5.9 x 10 8 


0237-233 


CS 1-0 


5.9 x 10 8 


0404+768 


NH3(1,1) 


1.2 x 10 11 


0428+205 


OHA1667 


5 x 10 9 




CO 1-0 


1.3 x 10 10 


0941-080 


CO 1-0 


2.0 x 10 10 


2352+495 


OHA1667 


2.8 x 10 9 



Note. - (1) B1950 IAU Source 
Name. (2) The observed transition. 
(3) The 3a upper limit on the column 
density in that line. (4) The equiva- 
lent upper limit to the molecular hy- 
drogen column density. (5) The es- 
timated upper limit to the molecular 
gas mass. For the limits on absorp- 
tion we assume a size of R = 3 kpc. 



